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^ ■ Abstract 

Based on realistic calculations for the nonmesonic decay rate of ^Hwe demon- 
strate, that in principle it is not possible to measure the total n- and p- induced 



decay rates and as a consequence T n /T p for that lightest hypernucleus. For 



the nonmesonic decay process the calculations are performed with modern YN 

forces based on various meson exchanges and taking the final state interaction 

among the three nucleons fully into account. Our findings have consequences 

also for the interpretation of experimental T n /Y p ratios for heavier hypernu- 

clei where severe discrepancies exist to theoretical r n /T p ratios. 
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I. INTRODUCTION 



There is a longstanding discrepancy between the theoretical ratio of the total neutron 
induced nonmesonic decay rate r„ to T p , the total decay rate for the proton-induced non- 
mesonic decay rate of various hypernuclei to experimental data Jl| . The experimental values 
are typically around 1 except for the very light hypernucleus ^He ]7|||, while theoretical 
evaluations lead to 0.05 - 0.2. The experimental value for T n is estimated either from neu- 
tron measurements and/or deduced from the measured values of the total nonmesonic decay 
rate T nm and of T p as 

r = r — r (i) 

L n — 1 n rn 1 p V / 

Apparently this relation can not be strictly true due to interferences. The quantity T p is 
determined experimentally from measuring single proton spectra and assuming that those 
protons are generated by the p-induced decay. Again this can not be strictly true since 
the n-induced decay leaves behind spectator proton(s) and final state interactions can carry 
momentum from neutrons to protons. We shall shed light in this article on those critical 
issues. On the theoretical side one faces the nuclear many body problem. Rigorous solutions 
based on realistic modern baryon-baryon forces are not in sight. Therefore shell model 
pictures supplemented by Jastrow type two-body correlations are typically being used and 
final state interactions are established by optical potentials. It appears difficult to estimate 
quantitatively the uncertainty of the theoretical predictions. In such a situation a view on 
very light systems is of increasing interest. In the 3-baryon system bound and scattering 
states can be rigorously gained based on modern realistic baryon-baryon forces p|. Therefore 
uncertainties about the quality of the hypernucleus wavefunction and final state interactions 
are absent. In the four-body system first rigorous solutions for bound states (^H and ^He) 
already appeared 0. The mesonic and nonmesonic decays of ^Hhave been calculated [pJQ 
but there are only few data to compare with. Some mesonic decay rates for ^Hfor which 
data are available agree rather well with that theory. Though there are state of the art 
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calculations no data are available for the very small nonmesonic decay rates of ^H. We 
would like to use in this article that theoretical insight to throw light on the questionable 
issues mentioned above. In [0] we found that the nonmesonic decays of |Hleading to a final 
deuteron and a neutron are suppressed by about a factor 10 with respect to the full breakup 
processes. Therefore we shall neglect those two-body fragmentation decay channels of ^Hin 
the following - except for pointing out that there a separation of n- and p- induced decays 
is clearly impossible. This is already evident from the fact that T^ +d + T p +d = 0.39xl0 7 s _1 , 
whereas the total n+d decay rate Y n+d = 0.66xl0 7 s -1 . Clearly there is a strong interference 
between the n- and p-induced decays. 

The exclusive differential n+n+p decay rate has the form 

* m,mi,m,2,m3 

M 2 k 2 k 2 
1^(2^ + ^-4)1 

Here k% and k% denote the directions of two detected nucleons (see || for further information). 
We have shown in that there are regions in phase-space which are populated by n- 
and p- induced decays and therefore an experimental separation for those contributions is 
impossible. But there are also regions in phase-space which are rather cleanly populated 
by either n- or p- induced processes, but not both. Therefore one has to be satisfied with 
certain fractions of T n and T p , defined by integrations over certain subregions of the total 
phase-space. In this manner one can measure n- and p- induced process separately. The 
fact that certain parts of the phase-space are populated by both processes coherently makes 
it obvious that by no means one will be able to access experimentally T n and T p separately. 
Nevertheless we would like to demonstrate this explicitely in the approach to T p which 
is being used for heavier hypernuclei JTJ. There one investigates the semiexclusive decay 
process in which only one proton is detected. Therefore we shall study in this article the 
single differential decay rate dT/dE p and in addition also dT/dE n and investigate whether 
they can be separated into n- and p- induced contributions and whether certain energy 
ranges are dominated by one or the other process. 
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Our results are based on rigorous solutions of the Faddeev equations for ^Hand the 3N 
final scattering states. We use the YN Nijmegen potential || which includes A-S conver- 
sion. It turned out that this potential produces the experimental ^Hbinding energy without 
further adjustment 0. For the NN forces we used the Nijmegen' 93 potential ||. We expect 
no dependence on the choice among the most modern NN potentials. For the hypertriton 
this has been verified. The importance of the final state interaction is demonstrated by also 
presenting results where the 3N scattering state in the nuclear matrixelement occurring in 
Eq. (^D is replaced by 3N plane wave states. This extreme approximation will, like in ||, be 
denoted by symmetrized plane wave impulse approximation (PWIAS), whereas the calcula- 
tion with final state interaction will be called "FULL". In Fig. 1 we show dT/dE n , dT n /dE n 
and dT p /dE n in PWIAS. The quantity dT / dE n has two peaks, one at very low neutron 
energies and one close to the maximal possible neutron energy. The peak at the higher 
energy is fed by the n- and p- induced processes as is obvious from the corresponding peaks 
in dT n /dE n and dT p /dE n . Clearly in both processes a high energetic neutron is produced. 
Surprisingly for us dT n /dE n + dT p /dE n sum up to dT/dE n with an error smaller than 5 %. 
The interference terms are therefore numerically very small. For very small neutron energies 
dT n /dE n dies out, since the n- induced process creates mostly high energetic neutrons. The 
p- induced process, however, dT p /dE n , exhibits a strong peak at very low neutron energies, 
which is caused by the (spectator) momentum distribution of the neutron in |H. Clearly a 
measurement of the decay rate dT/dE n as a function of the neutron energy will not allow 
to separate the n- and p- induced processes - except at very low neutron energies, where 
the energy distribution of the neutrons, however, is not determined by the A-decay process. 
That picture does not change qualitatively if one turns on the final state interaction as 
can be seen in Fig 2. Quantitatively, however, the rates are quite different. We can see a 
reduction factor of about 2 and the neglection of FSI would be disastrous in a quantitative 
analysis of data. Now the sum dT n /dE n + dT p /dE n equals dT/dE n only within about 12 
%. 

The situation for a separation of n- and p- induced processes appears somewhat more 
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favourable if one regards the single particle decay rates as a function of the proton energy. 
Our results are shown in Fig. 3 for PWIAS and Fig. 4 for the "FULL" calculation. For large 
proton energies nearly all protons result from the p- induced process : dT/dE p m dT p / dE p 
in case of PWIAS. The quantity dT n /dE p can not produce high energetic protons except due 
to FSI and this is indeed visible by comparing Figs. 3 and 4. dT n /dE p exhibits, however, 
the very low energetic proton peak from the spectator proton in |H. Also note again the 
reduction factor of about 2 caused by FSI. 

Let us now quantify the question, whether integrated proton distributions can provide 
a good estimate for Y p . Clearly the very low energetic peak should be excluded and one 
has to start integrating dT/dE p from the highest possible proton energy E™ ax , downwards. 
Thus we compare the integrals 



as functions of E p . The results are displayed in Figs. 5 and 6 for PWIAS and FULL. We see 
that in the case of PWIAS down to about E p ks 50MeV the two curves T(E P ) and Y P (E P ) are 
close to each other within less than 5% and only then start to deviate strongly. While T p (E p ) 
flattens out and approaches T P =T(E P = 0), T(E P ) receives contributions from the n-induced 
process. The situation is not so favourable, however, for the case FULL. Around E p =Q0 
MeV the relative deviation \T p (E p ) —T(E p )\/T p (E p ) is about 10 % and increase to about 20 
% around E" p =15MeV. Below that the deviation increases up to 30 %. Note also the relative 
factor of about 2 between PWIAS and FULL. We have to conclude that an estimate for 
T p from dT/dE p is only possible within an error of about 30 %. If one is satisfied with a 
fraction of T p the error can be reduced to about 10 %. 




(3) 




(4) 



and 




(5) 
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A well defined manner to receive information on the p- and n-induced decays separately 
is to use the differential decay rate of Eq.(Q) as described in 0. There are certain regions 
in phase space which are populated only by the p- induced decay and others which are 
populated only by the n-induced decay. In this manner one does not get the total T p or T n , 
but at least well defined fractions thereof. 

Now we would like to address the question, whether T n can be found via Eq. ([!]) in 
case of |H. This is a pure theoretical issue since, as we just demonstrated, T p can not be 
measured for |H. Surprisingly enough Eq. ([!]) is valid. As seen from Table V in we have 

T F n ULL = 0.17 x 10 8 , T FULL = 0.39 x 10 8 , T FULL + T FULL = 0.56 x 10 8 (6) 

That sum has to be compared with Y FULL = 0.57 xlO 8 , which treats the full process correctly 
as a coherent sum of the n- and p- induced decays. These numerical results validate Eq.(|l]) 
in the case of ^H. 

Finally we note that our theoretical result for the ratio of the total n- and p-induced 
decay rates in case of ^His T n /T p = 0.44. 

Since T p in the case of ^Hcan not be measured, it appears advisable to concentrate 
directly on dT/dE p and dT /dE n and compare those distributions to theory. This is an alter- 
native to the above mentioned exclusive processes. While measurements of the nonmesonic 
decay of |Happear to be far away, data for the four-body hypernuclei already exist f7|,|8| 
and theoretical predictions can be expected to come up in the near future. This will then 
allow interesting tests of the nonmesonic decay matrixelements, which will be based on real- 
istic four-body wavefunctions and various meson-exchange operators [@,|1(J > which drive the 
nonmesonic decay process. 
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FIG. 1. The single neutron decay rates dT/dE n , dT n /dE n and dT p /dE n in PWIAS as a 
function of the neutron energy E n . A separation in n- and p-induced processes is not possible. The 
peak at very low E n 's shows directly the momentum distribution of the neutron in ^H. 
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FIG. 2. The same as in Fig. 1 for the FULL calculation. The peak at very low E n 's is now 
also influenced by final state interactions. 
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FIG. 3. The single proton decay rates dT/dE p , dT n /dE p and dT p /dE p in PWIAS as a function 
of the proton energy E p . Now a separation in n- and p- induced processes would be possible for E p 
larger than about 50 MeV. The peak at very low E p s shows directly the momentum distribution 
of the proton in ^H. 




FIG. 4. The same as in Fig. 3 for the FULL calculation. The final state interaction causes 
now small contributions of high energetic protons resulting from the n-induced decay. Also the 
peak at very low E p s is now influenced by final state interactions. 
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FIG. 5. The integrated single proton decay rates according to Eqs.(H)-(||) for PWIAS. For 
E p >50MeV T P (E P ) &T(E P ). 




FIG. 6. The same as in Fig. 5 for FULL. Now the influence of the n- induced decay does not 
allow to estimate T p (E p ) by T(E P ). 
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